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Center clusters in full QCD at finite temperature and background magnetic field 
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We study the center structure of full dynamical QCD at finite temperatures and nonzero values 
of the background magnetic field using continuum extrapolated lattice data. We concentrate on 
two particular observables characterizing center clusters: their fractality and the probability for 
percolation. For temperatures below and around the transition region, the fractal dimension is 
found to be significantly smaller than three, leading to a vanishing mean free path inside the cluster 
structure. This finding might be relevant for center symmetry-based models of heavy-ion collisions. 
In addition, the percolation probability is employed to define the transition temperature and to map 
out the QCD phase diagram in the magnetic field-temperature plane. 


I. INTRODUCTION 

Quantum chromodynamics (QCD) is the theory de¬ 
scribing strongly interacting matter. QCD predicts the 
existence of a finite temperature transition that separates 
the low-energy confined regime and the deconfined quark- 
gluon-plasma (QGP) phase. The properties of this tran¬ 
sition are relevant for the evolution of the early universe 
and are also probed by contemporary heavy-ion collision 
experiments, both at RHIC and at the LHC. 

Following the conjecture that the deconfinement tran¬ 
sition in the gluonic sector is related to the magnetic 
transition of a corresponding spin system DU. and the 
finding that the latter can be understood in terms of 
cluster percolation [3], it was proposed that the gluonic 
field configurations of QCD can be characterized by cen¬ 
ter clusters and that the deconfinement transition may 
be understood as a percolation phenomenon [4]. In this 
description confinement manifests itself in small and un¬ 
correlated clusters, while the deconfined regime exhibits 
a large cluster that percolates and induces long-range 
correlations. The center structure of the QGP was also 
incorporated in models of heavy-ion collisions Di and 
was argued to explain various properties of the plasma 
phase including its low shear viscosity and high (color) 
opacity [6]. The main ingredient in this kind of models 
is the scattering of partons on the cluster walls, charac¬ 
terized by a mean free path. 

Besides the temperature, another parameter relevant 
for heavy-ion phenomenology is the background (elec¬ 
tromagnetic field generated by spectator particles in off- 
central collisions. Strong magnetic fields are also thought 
to have existed in the early stages of the universe and thus 
their effects on the QGP are of interest for cosmology as 
well. For recent reviews on the role of magnetic fields for 
strongly interacting matter see, e.g., Refs. Em¬ 
in this paper we perform numerical lattice simulations 
to study center clusters in 2 + 1-flavor QCD and deter¬ 
mine their response to nonzero temperatures and back¬ 
ground magnetic fields. We confirm that the clusters are 
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not three-dimensional objects but instead have a frac¬ 
tal nature, as has already been observed in pure gauge 
theory (see, e.g., Ref. 0). We demonstrate that as a 
consequence of this fractality the mean free path inside 
the clusters vanishes for temperatures and magnetic fields 
relevant for heavy-ion phenomenology. Furthermore, we 
propose a new observable for determining the transition 
temperature in full QCD and use it to map out the phase 
diagram in the magnetic field-temperature plane. 

II. CENTER CLUSTERS 

The concept of center clusters relies on the center sym¬ 
metry of pure gauge theory, formulated in Euclidean 
space-time at a nonzero temperature T. Center sym¬ 
metry denotes the invariance of the action under topo¬ 
logically non-trivial transformations g. These - unlike 
normal gauge transformations - are only periodic up to 
a constant twist, g(x , t+l/T) = zg(x, t ) in the Euclidean 
time-like direction m- Here, z belongs to the center 

Z 3 = {l,e _27ri/3 ,e 27ri/3 }, (1) 

of the gauge group SU(3). While the confined phase is 
center symmetric for pure gauge theory, in the decon¬ 
fined phase this symmetry is spontaneously broken. The 
corresponding order parameter is the expectation value 
of the Polyakov loop, defined on the lattice as 

p=yY, Tr ( 2 ) 

X t 

where the non-Abelian vector potential is represented 
by group elements C/ M = e xaA e , and V denotes the spatial 
volume of the system. For pure gauge theory the expec¬ 
tation value of P vanishes below the transition temper¬ 
ature T c and selects one of the center sectors ([T]) above 
the transition. In pure SU(3) gauge theory this decon¬ 
finement transition is of first order mm- 

The presence of dynamical quarks modifies this picture 
slightly: the fermion determinant breaks center symme¬ 
try explicitly and always favors the trivial center element 
1 (see, e.g., Ref. m)- However, this explicit breaking is 
rather mild and the Polyakov loop can still be used as an 
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approximate order parameter. The corresponding decon¬ 
finement transition is no real phase transition but merely 
an analytic crossover HUES]. For a pedagogical intro¬ 
duction to center symmetry and the Polyakov loop, see 
Ref. [TB] , 

Although the expectation value of P is - due to the 
explicit breaking - always real, it turns out that there 
are local domains in space, in which the Polyakov loop 
points towards one of the three center sectors EZH25]- 
The corresponding local Polyakov loops L(x) read 

L{x) = TrR U±(x, t), P=^L0r). (3) 

t X 

Below T c , all three sectors are (almost) equally rep¬ 
resented, giving rise to a cancellation and an (almost) 
vanishing average Polyakov loop P. This is visual¬ 
ized in Fig. |T| where the histogram of the local phase 
9{x) = arg L(x) is shown for a typical low-temperature 
configuration. For temperatures above T c , the real sec¬ 
tor 6 ss 0 becomes dominant (also included in Fig. |T]) 
and induces a large real average Polyakov loop. This pic¬ 
ture of center clusters has been studied in pure gauge 
theory with two mm, with three HEM and with 
four colors (21j . while preliminary results for dynamical 
quarks have been obtained in Refs. {221 G23]. (For visu¬ 
alisations of the clusters, see Refs. {HJ US]-) We men¬ 
tion that while the change in the distribution of arg L{x) 
is essential for the deconfinement transition, the mod¬ 
ulus |L(x)| was found to play no relevant role in this 
respect E® 

Besides the distinct population of the three sectors 
below and above T c , there is another pronounced dif¬ 
ference between the confined and deconfined regimes. 
While the clusters are small below T c , they percolate and 
span across the total volume above the transition region. 
In this sense the deconfinement transition becomes very 
similar to the percolation phenomenon in a three-state 
spin system. To give the center clusters a precise defi¬ 
nition that conforms to this picture, we need to impose 
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FIG. 1. Histogram of the local Polyakov loop phase below 
and above the transition temperature and the definition of 
sector numbers according to Eq. [4] 


a filter on the local phases 9{x) that discards sites lying 
far from center elements. Specifically, to each site x we 
assign a sector number n(x) £ {—1, 0,1} in the following 
manner m- 

( +i for e e [ § + <5 , 7T — <5], 
n(x) = { 0 for 0e[- f +8, J ■ /. 

[ — 1 for 6 £ [—7T -P 5 , — | — <5], 

(4) 


Here, / £ [0,1) is a free parameter, which removes “un¬ 
decided” sites, i.e., those that lie close to the minima of 
the distribution H(9), see Fig. [T] In the following we will 
refer to / as the cut parameter. The center clusters are 
then constructed in the following way: two neighboring 
sites x and y belong to the same cluster if their sector 
numbers are the same, that is, if n(x) = n(y). This di¬ 
vides space into domains where the local Polyakov loop 
points towards one of the three center elements. 

We emphasize that a nonzero cut parameter is neces¬ 
sary to interpret the deconfinement transition as a perco¬ 
lation phenomenon. Indeed, at / = 0, the center clusters 
would percolate already at low temperature^] By intro¬ 
ducing f ^ 0 and discarding sites lying far from center 
elements, the clusters are made thinner and percolation 
is delayed to set in only around T c . This way, the con¬ 
fined phase exhibits clusters with finite size, while in the 
deconfined phase there is one percolating cluster, as was 
demonstrated in pure gauge theory 0 Eni [20] . Note that 
similar thinning techniques (cf. Ref. 0 ) to reduce the 
cluster size are necessary in different contexts as well, 
e.g., for the magnetic transition in the Potts model m 
or for the droplet description of the Ising model Ea¬ 


rn. RESULTS 

The results presented below are based on the gauge 
configurations generated in Refs. [25H3T] at various values 
of the temperature, of the magnetic field B and of the 
lattice spacing a. These ensembles have been produced 
using the Symanzik tree-level improved gauge action and 
2 + 1 flavors of stout smeared rooted staggered quarks 
with physical masses. Details of the simulation setup and 
of the algorithm can be found in Refs. [28] [32J [33]. In the 
following we consider the stout smeared gauge links for 
calculating the local Polyakov loops. 

The vacuum configurations (corresponding to T ~ 0, 
B = 0) with several lattice spacings are used to set the 
cut parameter /(a) in a consistent manner. At finite 


1 To see this, note that in random percolation theory, the critical 
probability for a three-dimensional cubic lattice is p c rj 0.31 < 
1/3 m- Thus, even if the local Polyakov loops are completely 
random (i.e., the center sectors are equally populated) such that 
p = 1/3, each of the three sectors will percolate on average. For 
an explicit demonstration of this effect see Refs. mm ■ 
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temperatures we consider Nf x N t lattices and employ 
the fixed-TV* approach to vary the temperature. That 
is to say the temperature T = (7V t a) _1 is changed by 
tuning the lattice spacing a for a fixed lattice geometry. 
In this approach the continuum limit corresponds to the 
limit Nt —> oo at a given temperature. The magnetic 
field is chosen to point in the ^-direction and enters the 
simulation setup via its quantized flux, 

= eB ■ ( aN s f = 6nN b , N b e Z, (5) 

where the magnetic field is measured in units of the el¬ 
ementary charge e > 0. Due to flux quantization, an 
interpolation of the data at fixed N b is necessary to ob¬ 
tain results as a function of eB. For further details on 
the implementation of the magnetic field see Ref. [25]. 


A. Scale setting 



FIG. 2. The radius R of the largest cluster as 
a function of the cut parameter / for our zero- 
temperature ensembles with various lattice spac- 
ings. The dashed line indicates the prescribed clus¬ 
ter radius R = 0.49 fm. 


To set the cut parameter unambiguously additional 
physical input is necessary. A possible way to set / is 
to prescribe the value that the physical radius R of the 
largest cluster should take at low temperatures [9, 20] , 
For a cluster of size s, we define the radius R by the 
mean squared deviation of the sites r; in the cluster from 
its center of mass Rcm: 


Having fixed the precise definition of the clusters - 
i.e., the dependence of the cut parameter on the lattice 
spacing - at T = B = 0, we proceed to determine various 
properties of the clusters at nonzero temperatures and 
nonzero background magnetic fields. 

B. Fractality and the mean free path 


r 2 = 1 y>i - rcm) 2 . (6) 

To put this implicit prescription into practice we need 
to search for the value of / where the largest cluster 
has the desired radius. This procedure is visualized in 
Fig.[2]for various zero-temperature lattice ensembles with 
different lattice spacings a. Reading off the intersec¬ 
tion of the R(f) curves with the prescribed radius of 
R = 2.5 = 0.49 fm determines the scaling relation 

/(a). Note that for / —► 1, all sites are removed and, 
thus, the radius shrinks to zero, while for / = 0, the 
largest cluster fills the total volume so that R equals half 
the linear lattice size. The value R = 0.49 fm which we 
chose for setting / corresponds to a typical hadronic size 
relevant for the low-temperature confined regime. Note, 
however, that we are free to choose different radii as well. 
The subsequent analysis is performed using various val¬ 
ues 0.35 fm < R < 0.5 fm. 

The so obtained dependence /(a) is shown in Fig. [3] 
for various values of the fixed radius R. The curves all 
have positive slopes, as expected: for finer lattices the 
cluster radius in lattice units R/a has to be larger so 
that the radius in physical units R = a ■ R/a remains 
fixed. Thus, for smaller a the clusters must be made 
larger (in lattice units) via decreasing /. In the following, 
the interpolation of the /(a) curve will be used to set 
the cut parameter (for a few N t = 10 simulation points 
at high temperature, a controlled extrapolation is also 
necessary). 


We continue the analysis by demonstrating the fractal 
nature of the clusters. To this end, we employ the box¬ 
counting method to define the fractal dimension da ■ This 
approach is based on the scaling 

N(s) oc s -dn , (7) 


of the number N of boxes of linear size s necessary to 
cover a given cluster. This method was applied and com¬ 
pared to different definitions for pure gauge theory in 


Ref. 5 • 

Fig. 4] shows the fractal dimension of the largest clus¬ 
ter as a function of the temperature for several differ- 



a [fm] 

FIG. 3. The cut parameter as a function of the 
lattice spacing for various fixed cluster radii. 
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T [GeV] 


FIG. 4. The fractal dimension as a function of the 
temperature for three lattice spacings. The zero- 
temperature cluster radius is fixed to R = 0.49 fm. 


ent lattice spacings. At T = 113 MeV, where five lat¬ 
tice spacings are available, the a —> 0 extrapolation gives 
da = 1.9(1) in the continuum limit. For higher tempera¬ 
tures we find that three lattice spacings do not suffice for 
a controlled continuum extrapolation of this observable. 

In center cluster-based models of heavy-ion collisions a 
relevant parameter is the mean free path of partons inside 
the clusters. It is defined as the average distance that 
the parton can move without scattering on the cluster 
walls. To translate this notion into our setup, we consider 
the following procedure. For each site s inside a cluster, 
we count the number n l s of sites one can move in the 
direction i without reaching the boundary of the cluster. 
The (average) mean free path is then given by 

A /=3 E A / } ’ A / ) = ^E<-«’ (8) 

i=x,y,z s=l 

where S is the total number of sites available for the 
clusters. 

Above we have seen that the clusters are not three- 
dimensional objects but fractals. In the pure gauge the¬ 
ory setting it was pointed out already in Ref. [5] that 
as a consequence of this fractality the mean free path is 
not related to the linear cluster size but is much smaller 
than that. Using a continuum extrapolation based on 
three different lattice spacings, we show that A/ is consis¬ 
tent with zero for T < 300 MeV in the continuum limit, 
see Fig. [5] The systematic error of the continuum ex¬ 
trapolation is estimated by comparing fits with different 
forms for the T- and A^-dependence of A/. On a finite 
lattice, the fractal pattern is not resolved on distances 
smaller than the lattice spacing, thus the mean free path 
is bounded from below by a. Indeed, Fig. [5] reveals how 
the finite N t results for A/ approach zero via nonzero val¬ 
ues. Note that as the temperature is increased further at 
fixed lattice spacing a, and the largest cluster becomes 
three-dimensional, A f will approach half the linear lattice 
size (i.e., it will diverge in the infinite volume limit). 



T [GeV] 

FIG. 5. The mean free path as a function of the 
temperature, for three lattice spacings and a con¬ 
tinuum extrapolation. 


The background magnetic field breaks rotational sym¬ 
metry and thus might induce an anisotropy in the di¬ 
rectional mean free paths A^ A , defined in Eq. (jsj) . The 
effect of B on the Polyakov loop (and, thus, on center 
clusters) is indirect and occurs through virtual quark 
loops. In strong magnetic fields these virtual quarks oc¬ 
cupy Landau-levels: they are free to move parallel to the 
magnetic field but are localized perpendicular to it. This 
anisotropy is expected to propagate in the gluonic sec¬ 
tor and appear in the orientation of center clusters as 
well, implying A^ > A^ = A^\ Another argument 
supporting this hierarchy is based on the finding [34] 
that the magnetic field reduces the string tension in the 
parallel but increases it in the perpendicular direction. 
Indeed, a reduced string tension implies enhanced cor¬ 
relations between distant Polyakov loops and, thus, an 
increased mean free path in a given direction. Interest¬ 
ingly, in the asymptotically strong magnetic field limit of 
QCD [35] the parallel string tension even vanishes and 
local Polyakov loops are independent of 0 [3Tj. There¬ 
fore, in this limit center clusters become tubes in the z- 
direction but are expected to retain their fractal nature 
in the x — y plane. Nevertheless, our largest available 
magnetic field eB = 3.25 GeV 2 is still well below this 
asymptotic limit. 

To determine whether the predicted anisotropy is 
present in the center structure we calculated the direc¬ 
tional mean free paths at eB = 3.25 GeV 2 . In accordance 
with the above expectation we observe A), to exceed the 
perpendicular mean free paths, although only by a few 
percent. For lower magnetic fields 0 < eB < 0.7 GeV 2 
the effect is found to be smaller than our statistical er¬ 
rors. In this range the main effect of the magnetic field 
turned out to be described by a shift of the transition re¬ 
gion towards lower temperatures. We discuss this effect 
in more detail in the next section. 
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FIG. 6. The percolation probability as a function of the tem¬ 
perature for three different values of the magnetic field. The 
zero-temperature cluster radius is fixed to R = 0.4 fm. 


FIG. 7. The transition temperature, defined according to 
Eq. @, as a function of the magnetic field for three different 
lattice spacings. The zero-temperature cluster radius is fixed 
to R = 0.45 fm. 


C. The QCD phase diagram 

Due to the crossover nature of the deconfinement tran¬ 
sition in full QCD, the observables sensitive to the transi¬ 
tion exhibit no singular behavior but are instead smooth 
functions of the temperature. An implication of this is 
that the transition temperature is not uniquely defined: 
different definitions may result in different values for T c . 

The most straightforward definition involves the in¬ 
flection point of the average Polyakov loop. However, 
this turns out to be numerically difficult to locate due to 
the slow and gradual rise of P with the temperature, cf. 
Ref. [35]. There have been proposals to circumvent this 
issue by considering, e.g., ratios of Polyakov loop suscep¬ 
tibilities that take well-defined values both well above 
and well below T c , see Ref. m- 

Here we propose a new method to define T c using cen¬ 
ter clusters. In terms of the center structure, the most 
substantial difference between the confined/deconfined 
regimes is the absence/presence of percolating clusters^ 
(A cluster is defined to be percolating if it spans across 
the lattice in at least one spatial direction. Thus, such 
clusters become infinitely large in the infinite volume 
limit.) The simplest choice reflecting the abrupt change 
of gluonic configurations in this respect is the percolation 
probability p^ mmrnm- It is defined as the proba¬ 
bility of having a percolating cluster and is thus bounded 


2 Note that this direct realization of the Svetitsky-Yaffe conjec¬ 
ture becomes considerably more involved for SXJ(N) theories with 
N > 4. Unlike for SU(3) - where the Polyakov loop effective ac¬ 
tion is constructed exclusively via L(x) of Eq. - for SU(4), it 
involves the trace of gauge links in representations with different 
dimensions (4 and 6) |21M38lf39] . The 4-dimensional representa¬ 
tion alone was shown to be insufficient to describe the deconfine¬ 
ment transition via percolation, since the clusters were found to 
become too thin towards the continuum limit ED- (This is in line 
with the expectation based on random percolation theory, where 
the equally populated sectors have probability p = 1/4 < p c , cf. 
footnote 1.) Here we constrain the discussion to TV = 3, where 
such complications are absent. 


as 0 < poo <1. In Fig. [6] we plot Poo as measured on 
the 24 3 x 6 lattices, showing the expected rapid increase 
around T c . 

Taking into account the limiting values of Poo at low 
and at high temperatures, respectively, the most conve¬ 
nient choice for defining T c is through the implicit equa¬ 
tion 

Poo(T c ) = 0.5. (9) 

This definition will be employed below to map out the 
phase diagram for nonzero magnetic fields. 

To demonstrate the effect of magnetic field^ on poo, 
Fig-Ill also includes the percolation probability for a few 
nonzero values of B. Clearly, the magnetic field increases 
Poo for all temperatures and, as a result, reduces the 
transition temperature. This is consistent with previ¬ 
ous determinations of T C (B) using chiral quantities |28|. 
To quantify this effect, we employed the definition (|9| to 
determine T c for a range of magnetic fields using three 
lattice ensembles with N t = 6, 8 and 10. Fig. [7] shows 
the so obtained T C {B), revealing that the results for all 
three lattice spacings fall on top of each other. The tran¬ 
sition temperature is found to decrease by about 10% up 
to eB = 0.75 GeV 2 . 

Above, the cut parameter was set by fixing the zero- 
temperature cluster radius to R = 0.45 fm. Also it is of 
interest how the results change if R is varied. We have 
performed the same analysis for different values of R. 
Fig- i shows the continuum extrapolated transition tem¬ 
peratures based on our three lattice spacings for three 
values of the cluster radius, R = 0.40 fm, R = 0.45 fm, 
and R = 0.49 fm. The net effect of decreasing R is to 


3 We found that there is no anisotropy in the percolation proba¬ 
bilities, even for our strongest magnetic field. Instead, B only 
induces a weak anisotropy over shorter length scales, as revealed 
by the hierarchy in the directional mean free paths discussed in 
Sec.|MB| 
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eB [GeV 2 ] 

FIG. 8. Continuum extrapolated transition temperatures as 
a function of the magnetic field for different zero-temperature 
cluster radii R. (The last three points on the far right are not 
continuum extrapolated but were obtained on our N t = 16 
ensemble.) 


shift the transition temperature up. This is to be ex¬ 
pected: the smaller the low-temperature clusters are, the 
stronger ordering in the local Polyakov loops (i.e., the 
higher temperature) is necessary for percolation to set in. 
Notice that R affects T c because of the crossover nature 
of the transition, i.e., because the percolation probability 
depends smoothly on the temperature (even in the infi¬ 
nite volume limit). The gradual enhancement of p x (T) 
around T c becomes a real jump in pure gauge theory [9], 
where the transition is of first order. In the latter case, 
the clusters start to percolate suddenly, so that finite 
changes in the low-temperature cluster radius R are not 
expected to affect T c . Therefore, the change in T c due 
to varying R gives a measure for the width (strength) of 
the deconfinement transition. 

It has recently been shown that the QCD phase di¬ 
agram exhibits a critical endpoint for extremely strong 
magnetic fields ED, where the crossover turns into a first 
order transition (see also Ref. [ID]). Fig. [8] also shows 
T c at a very larg^j magnetic field for the three low- 
temperature cluster radii. A further decrease in T c by 
about 20% can be observed, again in agreement with pre¬ 
vious findings based on other observables |3 11 . Moreover, 
the difference between the T c curves for the different radii 
decreases by about 50% from eB = 0 to eB ss 3.25 GeV 2 . 
According to our reasoning above this shows that the 
transition becomes stronger as the magnetic field grows 
and the predicted critical point is approached. 


IV. CONCLUSIONS 

In this paper we have presented first continuum ex¬ 
trapolated results for various observables related to the 
center structure of full dynamical QCD. Center clusters 
were identified using a consistent thinning technique in¬ 
volving one parameter (the cut parameter /) that is fixed 
by prescribing the cluster radius R at low temperatures. 

Using this prescription, the fractal dimension of the 
center clusters was shown to be significantly smaller than 
three. We demonstrated that this leads to a vanishing 
mean free path in the cluster structure over the range 
of temperatures 110 MeV < T < 300 MeV. We found 
that the presence of magnetic fields eB < 3.25 GeV 2 
does not change this result qualitatively - even at our 
strongest magnetic field the anisotropy in the cluster ori¬ 
entation remains below a few percent. Thus, for a broad 
range of temperatures and magnetic fields that are rele¬ 
vant for heavy-ion collision phenomenology, the contin¬ 
uum extrapolated mean free path vanishes. This finding 
suggests a limited applicability for models that build on a 
finite mean free path for scattering processes in the QGP. 

Furthermore, we proposed a method to define T c in full 
QCD using the percolation probability and employed this 
definition to determine the phase diagram for nonzero 
background magnetic fields. The results unambiguously 
show a reduction of T c with increasing B , in good agree¬ 
ment with the results obtained using other QCD observ¬ 
ables m mi- In addition, the variation of T c when 
changing the zero-temperature cluster radius R was ar¬ 
gued to measure the width of the crossover transition. 
This quantity was found to gradually decrease as B grows 
and the predicted critical endpoint at extremely strong 
magnetic fields is approached. Altogether, our findings 
demonstrate that the deconfinement transition in full 
three-color QCD can be described as a percolation phe¬ 
nomenon. The analysis of further observables and the 
discussion of finite volume effects will be performed in a 
forthcoming study m■ Finally, we note that generaliza¬ 
tions of the percolation picture to other gauge groups, 
e.g., SU (N) with N > 3, are non-trivial, and that more 
extensive research is required to adress their viability. 
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This is still well below the estimated critical magnetic field 
eB CEP = 10(2) GeV 2 [31]. 
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